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1. INTRODUCTION 

Power electronics play major roles in wind energy. They are responsible for several functions such 
as starting the system, regulating turbine speed, interfacing turbine with the grid, regulating real and reactive 
powers, and controlling the system during disturbances [1]. The doubly fed induction generator (DFIG) has 
many benefits. This is possible due to its ability to supply power at constant voltage and frequency while the 
rotor speed varies [2], [3], which allows better wind capture, compared to fixed speed wind turbines. 

The DFIG is a wound rotor induction generator, the stator winding of which are directly connected 
to the grid. The rotor winding is connected to the grid by means of a reversible static power frequency 
converter (inverter/rectifier) [4]. The high sensitivity of the DFIG to electrical faults reveals numerous 
difficulties in terms of conformance to specifications imposed by the electrical grid managers [2]. Defects 
appearing in an electrical system have various causes. They may be due to the aging of the system 
components (electronic components, machine coils) or to conditions of use or simply to a manufacturing 
defect whose effect would be imperceptible during commissioning [5], [6]. 

In a static converter feeding DFIG (two level inverter), a severe defect in an element of the energy 
conversion chain alone can cause the total or partial loss of the system control (DFIG+inverter). In this paper, 
Lissajous curves will be presented as a diagnosis method for detecting the inverter feeding doubly fed 
induction generator in the case of single and multiple short circuit switches faults. The proposed method has 
shown its effectiveness in several fields like the rotating machinery [7], faults diagnosis [8]-[10], it was also 
studied in the case of multiple open switches faults in three phase voltage inverter feeding DFIG presented in 
a previous paper [11]. 
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2. WIND TURBINE MODELING 
According to the theory of momentum and Bernoulli’s theorem, incident theoretical power due to 
wind Pwina is given by the equation as shown in (1). 


1. 
Pywind = 5 P *A* v? (1) 


Where A: Sweep area swept by the blade, 
P: Density of air, 1.225 kg/m, 
v: Wind speed, 
The expression of the mechanical power captured by the blade [12] is (2). 


Prurbine = = Cp g p FA # v? (2) 
Where Cp is the coefficient of power. 
According to [13], [14] this last coefficient is approximated by (3). 


q(à-—3) 
15—0.3ß 


Cy (A, B) = 0.4233 sin ( 


) — 0.0184(A — 3)B (3) 


Where À: specific speed, 
p: angle of the blades orientation. 
Then the aerodynamic torque is written (4) [15]. 


2 3 
C = Pturbine = 1 Cp *p*n*R" *V (4) 
aero Q 2 Q 


The multiplier assures adaptation of the slow speed of the blades to the fast speed of the machine by the gain 
G. The mechanical torque Cm and the speed of the machine Q, are presented respectively in (5) and (6): 


Cm = Caero/ GT (5) 
ipren 6) 


The fundamental equation of dynamics is given by (7). 


J a = as = Cm — Com — Cris = Cm — Com — f * Oy (7) 
Where J : Total inertia that appears on the generator shaft 
P : number of pole pairs 
0, : Electrical angular speed of the rotor 
Com : Electromagnetic torque 
Cvis : Viscous friction torque 
f : Coefficient of viscous friction 


3. MODELING OF DOUBLY FED INDUCTION GENERATOR 
Model of the asynchronous machine equations for the stator % and the rotor V, voltages in Park 
reference frame (d,q) are given by (8) [16]. 


d 
Vsa = Rs * lsa + qe Psa = Ws * Psq 


d 
Vaq = R; * Isq F Te (Psa) — Ws * Psa 


d 
Vra = Rr * Tra + dt (Pra) — (ws-wr) * Prq 


a (8) 
Vig =R, * Irq T AC — (ws-wr) * Pra 
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The electromagnetic torque equation is presented as shown in (9). 


3 Lm 
Cem = a i I (pra * Irq — Prq * Isa) (9) 


4. FIELD ORIENTED CONTROL 
In order to perform the DFIG active and reactive power by controlling dynamically and separately 
the flux and the torque, we choose orienting the stator flux along the direct axis and obtain: 


fee = Ps 
Pqs = 0 
M . 
Com = -P.T Pas: lar (10) 


Active and reactive powers can be written is being as (10). 


Fis. aD 


5. VOLTAGE SOURCE INVERTER SYSTEM 

Voltage source inverter (VSI) is an electronic device that relates direct current (DC) to alternating 
current (AC). VSI has benefits such as higher efficiency, reduced installation time, elimination of 
interconnect power wiring costs and reduced building floor space [17]. In a wind energy conversion system 
using a doubly fed induction generator, as can be seen from Figure 1 [18], a fixed DC input voltage is 
supplied to the inverter and a controlled AC output voltage is achieved by setting the on and off periods of 
the inverter switches. 

In this application, we are interested in the hypo-synchronous operation of the generator, so as to 
treat only the transfer of rotor power in one direction (grid to rotor). The grid side converter will therefore be 
a simple unidirectional rectifier. The DC bus voltage is therefore kept fixed and we are interested only in the 
control of the rotor side converter [11]. 


Figure 1. Wind energy conversion system using a doubly fed induction generator 


6. LISSAJOUS CURVES FOR DFIG DEFECTS DIAGNOSIS 

Lissajous curve is the trajectory of a point whose rectangular components have a sinusoidal 
movement, the x signal and the y signal are both described separately by different sine waves. So, by 
expressing y as a function of x, different patterns are created. In order to apply these curves to the DFIG rotor 
currents, it is necessary to use the Concordia transformation [19] on the three-phase rotor basic currents lape 
to obtain the rotor currents Concordia components (Iæ, Ig) such as: 


-2 a p" - 12 
ce B v3 | | ir 
2 ~ Z| Ure 
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7. SIMULATION RESULTS 

Simulation of the energy conversion chain was performed using MATLAB‘®- Simulink*® 
software. It allowed us to create multiple fault scenarios that could appear in the converter switches (single 
fault, double fault in the same leg or in different legs, possibility of cascades of defects of the same nature or 
different). 


7.1. Free-defects operation 

The Lissajous curve of doubly fed induction generator rotor currents Concordia components, as 
shown in Figure 2, has a shape of two circles of which the one inside has six knots: the number of the 
inverter switches. Subsequently, cases of defected inverter feeding the DFIG will be compared with this case 
reference. 
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Figure 2. Lissajous curve of a DFIG rotor currents during free-defects operation 


7.2. Faulty operation 

Various failures can impact the inverter during its operation, especially in terms of its so-called 
power components because of their fragility [20]. Two types of faults can be signalted: short circuit and open 
circuit faults [21], [22]. Note that a fault in an inverter feeding an asynchronous machine can cause 
significant consequences on the mechanical level: the ripple of the torque and the speed bring undesirable 
vibrations and annoying sound noises [5]. In some cases, the system is able to operate but with degraded 
performance (operates with imbalance in electrical and mechanical quantities and with low gain). In this 
paper, the rotor side converter faulty operations are considered in the case of multiple short-circuits switches 
faults that occur at the same leg and at different legs as can be seen from Figure 3 [23]. The short circuit fault 
condition occurs when one of the switching cell components (insulated gate bipolar transistor (IGBT) in our 
case study) stays closed constantly (S;=1). 


GRID INVERTER 


POWER gPs 


Figure 3. The six switches {5, to Sg} of the rotor side converter concerned by the faulty operation study 


Short-circuit faults can be considered as the most serious faults attacking the IGBT switches [24]. 
The current in the presence of this fault attains limits that can lead to the fusion of its chip or its connection. 
If the fault detection does not took place quickly (less than 10 microseconds), then the active IGBT switch on 
the same leg goes through the same phenomenon and therefore the entire leg of the inverter is shorted [25]. 
Note that this type of fault causes disturbed and unstable operation of the system due to the operating mode 
change during the application of the fault. 
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Indeed, short circuit fault can introduce a decrease in speed to a value almost zero in certain 
intervals with the presence of undulations, and increasing torque and currents to intolerable levels in other 
intervals. This allows us to say that this defect alone can cause enormous damage to the entire wind turbine 
conversion chain. Simulation results of the DFIG fed by a three-phase voltage inverter controlled in sine- 
triangular PWM are obtained. The Lissajous curves of rotor currents Concordia components indicate the 
inverter faulty situation. 


7.3. Inverter single switch fault 

The Figures 4 (a) and (b) show the rotor currents Lissajous curve in the case of different level single 
switch fault. When the fault occurs, the resulting shape is close to a circle with five knots with more 
concentration on one of the existing knots whose position depends on the defected switch location. In the 
case of one upper switch fault, the concentration is on the knot at the top Figure 4 (a). In the case of one 
lower switch fault, the concentration is on the knot at the bottom Figure 4 (b). 
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Figure 4. Lissajous curves of rotor currents Concordia components (a) case of a DFIG fed by S4 switch 
defected inverter and (b) case of a DFIG fed by S, switch defected inverter 


7.4. Inverter double switches fault in the same leg 

The Figures 5 (a)-(c) show the rotor currents Lissajous curve in the case of the same leg double 
switches faults. When the fault occurs, the resulting shape is close to a circle with four knots forming a 
butterfly with knot inside Figure 5 (a) until Figure 5 (c). The obtained shape changes its orientation according 
to the defected switches location. 
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Figure 5. Lissajous curves of rotor currents (a) case of a DFIG fed by S4 and S, switches defected inverter, 
(b) case of a DFIG fed by S3 and S, switches defected inverter, and (c) case of a DFIG fed by S; and S, 
switches defected inverter 


7.5. Inverter double switches fault in different legs 
The Figures 6 (a)-(c) show the rotor currents Lissajous curves in the case of different leg double 
switches fault. When the faults occur, the resulting shape is close to a circle with four knots instead of six 
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Figure 6. The circle missed knots changes its orientation according to the defected switches location. In the 
case of two upper switches fault, the missed knots appear at the bottom Figure 6 (a). In the case of two lower 
switches fault, the missed knots appear at the top Figure 6 (b). In the case of one upper switch and one lower 
switch faults, the two missed knots appear one at the top and the second at the bottom as shown in Figure 6 (c). 
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Figure 6. Lissajous curve of rotor currents (a) case of a DFIG fed by S, and S} switches defected inverter, 
(b) case of a DFIG fed by S, and S, switches defected inverter, and (c) case of a DFIG fed by S4 
and S,switches defected inverter 


8. CONCLUSION 

The work presented in this paper is devoted to the analysis, modelling and simulation of a variable 
speed wind turbine using a doubly fed induction generator DFIG. In the first part, we gave a brief description 
of the VSI, its operating principle and its advantages. In the second part, we studied simulation results of the 
inverter feeding the DFIG in the case of single and multiple short circuit switches faults. Several fault 
scenarios are applied and visualized using the method of Lissajous curves of DFIG rotor currents Concordia 
components which shows the defect exact location. In the case of a free-defects operation, the Lissajous 
curves have a shape of two circles of which the one inside has six knots (corresponding to the number of the 
inverter switches). The width of the circle connecting the 6 knots positions differs, the number of knots 
decreases, and the orientation of the circle changes depending on the number of the inverter defected 
switches and their position. The next work will deal with other powerful diagnostic methods that will be 
associated to the Lissajous curves method to determine the severity of the defect encountered and thus 
remedy this problem to avoid unwanted consequences. 
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